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ABSTRACT 

The total masses ejected during classical nova eruptions are needed to answer two 
questions with broad astrophysical implications: Can accreting white dwarfs be pushed 
towards the Chandrasekhar mass limit to yield type la supernovae? Are Ultra-luminous 
red variables a new kind of astrophysical phenomenon, or merely extreme classical no- 
vae? We review the methods used to determine nova ejecta masses. Except for the 
unique case of BT Mon (nova 1939), all nova ejecta mass determinations depend on 
untested assumptions and multi-parameter modeling. The remarkably simple assump- 
tion of equipartition between kinetic and radiated energy (E\^ n and E Ta( ±, respectively) 
in nova ejecta has been invoked as a way around this conundrum for the ultra- luminous 
red variable in M31. The deduced mass is far larger than that produced by any classical 
nova model. Our nova eruption simulations show that radiation and kinetic energy in 
nova ejecta are very far from being in energy equipartition, with variations of four orders 
of magnitude in the ratio E^ m /E va ^ being commonplace. The assumption of equiparti- 
tion must not be used to deduce nova ejecta masses; any such "determinations" can be 
overestimates by a factor of up to 10,000. We data-mined our extensive series of nova 
simulations to search for correlations that could yield nova ejecta masses. Remarkably, 
the mass ejected during a nova eruption is dependent only on (and is directly propor- 
tional to) -E r ad • If w e measure the distance to an erupting nova and its bolometric light 
curve then E ra & and hence the mass ejected can be directly measured. 

Subject headings: accretion, accretion disks — binaries: close — novae, cataclysmic 
variables — white dwarfs 
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1. Motivation 
1.1. SNIa 

Which type of star (or stars) is the progenitor of type la supernovae (SNIa)? The discovery, 
via SNIa, of the likely acceleration of the expansion of the universe inexorably drives us to the 
concept of dark energy. If correct, this breakthrough in our understanding of cosmology is pro- 
found... but significant uncertainties remain. One of the most perplexing of these uncertainties is 
that we still can't identify the progenitors of type la supernovae. Systematic changes in these pro- 
genitors over a Hubble time could systematically affect SNIa light curves in ways we can't predict; 
a worrying prospect for the brightest standard candles known. There is a broad consensus that if a 
white dwarf (WD) can somehow be made to accrete enough matter to exceed the Chandrasekhar 
mass, it will fuse carbon and explode, releasing ~ 10 51 ergs in photons and radioactive nickel-56 
(observationally-constrained quantities). The problem, of course, is that there are multiple possible 
donors able to dump matter onto a white dwarf in a binary system: brown dwarfs, white dwarfs, 
main sequence stars, red giants and supergiants. The absence of hydrogen in SNIa spectra is an 
important constraint, but unfortunately it doesn't definitively rule out any of the donor suspects. 
Just a tiny amount of accreted matter (~ 3 x 10~ 10 M© yr -1 )) - with or without hydrogen - can 
push a nearly-Chandrasekhar mass WD "over the edge" . Thus none of the hydrogen-rich donors 
noted above are excluded, at least in principle, if the accreting WD gains mass in secular fashion. 

WDs accreting hydrogen-rich matter usually undergo classical nova (CN) eruptions, which 
periodically eject mass. Knowing whether such WDs increase or decrease in mass as a result of 
these competing processes is equivalent to knowing if WDs in mass-transferring binary systems, 
with hydrogen-rich secondaries, can be the progenitors of SNIa. An effective way to determine if 
the mass of an accreting WD is secularly increasing would be to measure its rate of accretion and 
time between eruptions. The accreted mass must then be compared to the mass ejected (m e j) to 
determine if there is a net gain or loss in mass. Recent attempts to carry out this experiment with 
the recurrent nova ( RN) T Pyx, and the diffi culties with directly measuring m e j, are reported in 



Schaefer et al.1 (j2009T ) and lSelvelli et al.1 (|2008l l 



T Pyx is the best studied and prototypical RN. RN might give rise to SNIa; t he importance 



of kn owing whether T Pyx will (or will not) become an SNIa cannot be overstated. ISchaefer et al 



(|2009l ) have mustered strong evidence that T Pyx must have under gone a classical nova eruption 



around 1866, after a 2.6 Myr hibernation episode (IShara et all 1 19861 ) followed by a 750 Kyr phase 
of accretion at a rate of ~ 4 x 10 -11 M & yr _1 . Allowing for surrounding interstellar matter swept 
up by the expanding nova envelope, the 1866 classical nova ejecta mass should now total at least 
750 Kyrx ~ 4 x 10~ n M & yr _1 ~ 3 x 10 _5 M Q . The mass accreted and ejected since 1866, in 
five recorded recurrent nova outbursts, must be a few orders of magnitude less. If this scenario 
is correct then T Pyx is unlikely to become a SNIa, because the WD mass is secularly decreasing 
due to classical nova eruptions. A direct test of this scenario - and whether the prototypical 
recurrent nova will become an SNIa - would be an unambiguous measurement of the mass in the 
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ejecta surrounding T Pyx. A direct measurement of the total mass surrounding T Pyx is very 
challenging, as successive generations of ejecta catch up with and shock gas already in pl ace. Some 
shocked blobs appear and disappear on a timescale as short as a year (jShara et allll986l ). 



1.2. Ultra- luminous Red Novae 



Over the past decade it has been suggested that a few rare, but extra ordinarily lum inous 
eruptive variables are pr ototypes of a new class of astrophysical p henomenon rtHarwitl 119811). The 
"Red Variable" in M31 dRich et alJll98sh and (rMould et al.l[l99ol), V838 Mon (jBond et alJl2003h 
and the "Optical Transient" in M85 (jKulkarni et al.ll2007h and (|Rau et all 120071 ) reached absolute 
magnitudes as luminous as —10 to —12. This is brighter than classical novae, but fainter than 
supernovae. 



Classical novae typically eject ~ 1 — 10 x 1O _5 M0 during an outburst ( Yaron et al.ll2005l ). 
It is evident that mass ejections at least 100 times larger must have accompanied the outbursts 
of the Ultra-luminous Red Novae, as the ejecta remained optically thick to much larger distances 
from the sites of the explosions. As a result, the ejecta cooled to temperatures as low as 700 
Kelvins. Explaining these observations is challenging. Some of the a stronomical community sup- 
ports the hypothesis of "mergebursts" (the merger of a binary star) (jSoker fc Tvlendal 120031 ) and 
(jTvlenda fc Sokei 2006). The m ass ejected in such an event could easily exceed ~ 10~ 2 M Q . A com- 
peting model (jShara et gill 12010 ) posits that extreme classical novae (with very low (0.5M©)) WD 
masses and very high accreted envelopes (~ 1O _3 M0) can explain the Ultra-luminous red variables 
without invoking a new type of astrophysical phenomenon. An important numerical prediction of 
the extreme nova scenario is that the mass ejected in such an event is at most ~ 2 — 3 x 1O _3 M0. 
Thus a direct measurement of m e j in excess of ~ 10~ 2 M Q could eliminate one of the competing 
explanations for Ultra-luminous red variables. 



2. Measuring Ejecta Masses 



The preceding section emphasizes how desirable it would be to accurately measure the masses 
of ejecta from classical novae. The ejected masses of classical novae have long been debated and esti- 
mates have widely varied, depending on the methods employed to interpret observations. Estimates 
of the density of ejected gas (from the presence or absence of forbidden emission lines) and size of 
the ejecta (from the angular size and velocity of the ejecta) lead to the simplest possible estimated 
ejecta masses. However, the oft-mad e simple assumpt io ns of homogeneous, spherically sy mmetric 
mass ejection are demonstrably false (jShara et allll997l ). (jWade. Harlow h Ciardulldl2000i ). Filling 
factors very different from unity are probably c ommon, as are dense clumps in outflowing nova 
winds (j Williams fe Shafterll2004l ). iFerlandl (|1998l ) has emphasized that nova ejecta masses are sys- 
tematically underestimated because gas emitting in any particular region of the electromagnetic 
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spectrum (say the ultraviolet) can have very different properties from gas emitting in another region 
(say the thermal infrared). A region of gas may be an efficient emitter in one regime and totally 
invisible in another. These factors can easily lead to uncertainties of orders of magnitude in the 
estimated ejected masses. 

Could any or all of these factors - dumpiness, variable filling factors, variable ionization, and 
especially non-spherical ejection, - render all of the detailed 1-D nova models too doubtful to 
trust? While highly non-spherical ejection might compromise the ID models, the other factors are 
challenges in interpreting the observations rather than the models. Fortunately there is one very 
dependable measurement of a nova's ejected mass, whic h serves as a critical "sanity check" for 
extensive simulations of nova eruptions (jYaron et al.ll2005l ). ISchaefer &: Pattersonl (|1983l ) measured 
the period change of the eclipsing classical nova BT Mon (nova 1939) as a result of mass loss 
from the binary system. The modern period of this binary is 40 parts per million longer than 
before 1939. Reasonable assumptions concerning the angular momentum carried off by the ejected 
material provide a dynamical measure of the ejecta mass. The mass of the WD in BT Mon has 
been measured with rather high precision: 1.04 ± 0.06M Q (jSmith et allll998l ). 



Repeating ISchaefer &; Pattersonl (ll983l Vs derivation but with the more precise white dwarf 

4 x 1O- 5 M , to 



Smith et all (Il998h yields the BT Mon ejected 



mass: 



mass determinations of 

within a factor of two or so. A heroic reconstruction of t he light curve of BT Mon from Harvard 
archival photographic plates (jSchaefer &: Pattersonl Il983l ) suggests a time to decline, by two (£2) 
and three (£3) magnitudes, of 140 and 190 days, respectively. If correct, these would be useful 
constraints for testing nova simulations. Unfortunately, the maximum light of BT Mon w as almost 
certa in ly missed (wh en the star was in the da ytime sky), based on the outburst spectra (jWhipple 



1940j), (jSanfordlll940l ). (jSwings h Struvdll94ll ). These spectra demonstrate that the outburst must 
have occurred when the nova was at least 5.4 magnitudes brighter than when the photographic 
record began. Thus the t<i and £3 above do not apply to maximum light; the correct values must 
be considerably smaller, but we can't measure them directly. Fortunately though, the full width s 
of the emis sion lines were , well observed, c orresponding to velocit ies of 1500 km/s (I Whipple! 1 1 9401 ) . 



2100 km/s (|Sanfordlll94ol ) and 1730 km/s ([Swings k. Struvil94lh . We adopt the maximum ejected 
velocity as half the average of these three observations, viz. 900 km/s. 



Each of the nova simulations of lYaron et al.l (|2005l ) predicts an ejected mass and a maximum 
ejection velocity for a given white dwarf mass, mass accretion rate and luminosity. These quantities 
vary quite smoothly between simulations, so it is easy to interpolate. It is reassuring that, for a cold 
WD (10 Million Kelvins core te mperature) of LOMq , accreting at a rate of ~ 3 x 10~ 10 M@ yr~ , 
the 1-D nova simulation grid of lYaron et al.l ( 20051 ) predicts an ejected mass of ~ 7 x 1O~ 5 M 
and a maximum ejection velocity v of 1000 km/s, in very good agreement with the observations. 
(A further check would be possible if the WD luminosity were directly measureable, but because 
of ongoing accretion this isn't possible. There is also no easy way to check on the BT Mon WD 
core temperature, but 10 Million Kelvins is certainly plausible). Unfortunately not a single other 
dynamical mass for nova ejecta has ever been measured, so the ejected masses of Ultra-luminous 
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red variables, recurrent novae and of all other novae remain uncertain. 



This unc ertainty is overcom e if one assumes equipartition of energy between radiation and 
motion (e.g., iMould et al.l (jl990l )) in nova ejecta. As the total radiated energy E ia ^ is readily 
obtained by integrating the light curve, and an average expansion velocity v av may also be obtained 
from spectrographic observations, the ejected mass is derived: m e j ~ 2E Ia d/v% v . If the assumption 
of equipartition of energy could be trusted then a simple and powerful tool would be available for 
testing nova theory, and ideas about SNIa progenitors. 



3. Non-Equipartition Between Kinetic and Radiated Energy 



We decided to test this ad-hoc and very important assumption via our extensive evolutionary 
calculations of classical novae. The hydrod ynamic, 1-D Lagrangian stellar evolution code used in 
all our studies is described in some detail (jPrialnik fc Kovetall995l ). It includes OPAL opacities, 
an extended nuclear reactions network comprised of 40 heavy element isotopes, and a mass-loss 
algorithm that applies a steady, optically thick supersonic wind solution (following the phase of 
rapid expansion). In addition, diffusion is computed for all elements, accretional heating is taken 
into account and convective fluxes are calculated according to the mixing length theory. The code 
finely subdivides the outer layers so that the radius at which optical depth becomes unity is well 
determined, as is the effective temperature. This, in turn, allows the calculation of the radiated 
energy at each time-step. 

Initial models were prepared for a range of WD-mass values and three temperatures by cooling 
WD models from higher temperatures. Each nova model was followed through several consecutive 
outburst cycles in order to eliminate the effect of the initial configuration. One typical cycle was 
then chosen as representative of each parameter combination. The kinetic energy E\^ m in each 
simulation was obtained by summing ^M c j5tv 2 , where M e j is the mass ejection rate, over the entire 
mass loss episode. Similarly the total radiated energy E IS ^ was computed by summing the radiation 
emitted at all time-steps from the model's photosphere, L5t. 

Results for 58 simulations, u sing the full range of no va white dwarf ma sses, temperatures and 
accretion, including the models of lYaron et al.l (|2005l ) and lShara et all (|2010l ) are shown in Figure [TJ 
£ , ra d/-E'kin values of 10 to 1000 are the norm. But we see that the ratio of radiated to kinetic energy 
can be as small as 1 or large as 10,000! Co nfirming this , a single, very detailed simulation of a 
solar mass white dwarf yielded a ratio of 10 (|Kovetall999l ). Figure [1] demonstrates that E ia ^/E^ n 
tends to decrease with increasing WD mass, but we need more information - the underlying WD 
luminosity and accretion rate - to uniquely determine the value of £7 r ad/ -^kin- If we could somehow 
measure the WD mass, its intrinsic luminosity, and the accretion rate before the eruption we could 
then use Figured] to uniquely determine -Brad/^kin- This information is simply not available. If it 
were, we could in fact determine the ejected mass directly from the grid of models. 



Mass estimates that simply assume that i^rad/^kin = 



1 will be in error by up to a factor of 
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10,000. Thus, assuming equipartition results in gross overestimates of the ejected masses derived 
from observations of ejection velocities and nova luminosities. We cannot emphasize too strongly 
that the assumption of equipartition between luminosity and motion in nova ejecta is unjustified 
and leads to extremely overestimated masses. This may help settle the controversy between the low 
ejecta masses advocated by nova modelers and the (sometimes) significantly higher masses claimed 
by observers. In particular, equipartition-based suggestions of 1O _2 M0 or even KHAf© for the 
ejecta of objects like M31-RV should be discounted. 

Even though the assumption of equipartition fails in predicting m e j, we are not forced to end 
the discussion on a negative note. In the following section we demonstrate that nova ejecta masses 
can be determined with readily available observational data. 



4. Nova Ejecta Masses from Ejection Velocities and E TS ^/E) 



kin 



In principle, if £7 r ad 7^ ^kim the nova ejecta mass may still be obtained from 

m ej = [E Tajd /(0.5v 2 )][E kin /E Tajd \, 



(1) 



where the first term is obtained purely from observations and the second term is provided by 
models. 

The (distance-independent) quantity v requires a series of spectra of a nova to determine v. 
Measuring E rad itself can be more challenging: an accurate distance to a nova and its bolometric 
light curve are requi red. T h e dist ance can be determined fr o m th e absolute magnitude - decline 
time relationships of ISharal (|198ll ) and iDownes Duerbeckl (|2000i ) , or the expansion parallax of 
the nova shell, or the distance of a host star cluster or galaxy in which a nova is located. In any 
such derivation it is essential to include the very large contributions o f ultraviolet and infrare d 
radiation from novae, typical ly seen after the pea k in visual brightne ss ([Gallagher &: Code! Il974l ) , 
([Gallagher fc Starrfieldlll97^ ). (|Geisel et al.lll97oh . ifaehrz et alJll98cl ). 



The pro blem is how to connect m odels with observations for as few observables as possible. 
As shown bv lPrialnik &: Kovetzl (|1995l ). nova outbursts constitute a 3-parameter family of events. 
Three independent observable parameters that span the parameter space may be taken to be the 
helium abundance Y and the metallicity Z in a nova's ejecta, and the t ime of decli ne of the 
nova luminosity £3 (or some other related observable like ejection velocity (|Sharalll98ll )), These 
parameters, the composition-related ones in particular, are not always easy to determine. Do Y or 
Z in a nova's ejecta control the value of E^/E^l In Figure[2]we have plotted -Erad/^kin versus Y 
and Z for 58 nova models. The resulting scatter-plot tells us clearly that neither Y nor Z correlate 
with EVad/i^km- The correlation with £3 is slightly better but not at all convincing. 

Is there some other observable parameter that correlates with E Iaj ^/ 'E^l As radiation pressure 
is the main driver of mass loss during a nova eruption, we might expect the mass loss rate M e j 
to behave similarly to Reimers's mass loss formula that is used in conjunction with red giant 
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winds: M e j oc LM/R. Since the photospheric radius and WD mass change little during the mass 
loss episode of a nova, it follows that M e j might be roughly proportional to the luminosity. If 
this were the case, then the total ejected mass m e j should be proportional to the product of a 
nova's luminosity and duration of mass ejection ~ £? ra d- Combined with E^ = 0.5m e jV 2 , these 
relationships suggest that i£ ra d/ -^kin should be proportional to v~ 2 . 

Our nova simulations can be used to test this heuristic suggestion. Figure O using the same 58 
simulations of Figure HJ shows the above suggestion to be correct, with remarkably small scatter, 
when we plot £J ra d/-Ekin versus the average ejection velocity v: 

log(£ rad /£kin) = -2.01og(«) + 7.6 (2) 

where the average ejection velocity v is measured in km/s . Equations (1) and (2) immediately 
lead to 

m ej = 6 x KT 18 £ rad (3) 

where the ejected mass is measured in gm and the radiated energy in ergs. Figure 4 is a plot of 
the 58 models' computed ejected mass versus radiated energy. The line through the data points in 
Figure 4 is equation 3, the key result of this paper. 

We thus derive the remarkable result that the mass of a nova's ejecta can be determined to 
within a factor of two (as seen from Figure 4) if we can measure just the time-integrated bolometric 
luminosity; spectra to determine v are not needed. It is still a significant task to determine both the 
distance and especially the UV through IR flux of a nova, but this information alone is sufficient 
to determine the long-sought masses of nova ejecta. The coefficient in equation 3 depends on our 
grid of nova models; its precise value should eventually be directly measured in a few novae with 
well-determined ejecta masses and total radiated energies. 

Unfortunately, BT Mon and its dynamically determined shell mass cannot be used to test this 
methodology. The nova reached peak brightn ess when it wa s in the daytime sky, it had faded at 



least five magnitudes before it was discovered (j Whipple! 1 194Q ). which was decades before IR or UV 



measurements were possible. Thus the total radiated energy (including the contributions to E ra ^ 
from the emitted UV and IR radiation) are unknown. Still, we point out that there is reason to be 



optimistic for the future. ISchaefer Pattersonl (|1983l ) succeeded in measuring the BT Mon period 



change because of the treasure trove of plates in the Harvard plate stacks. It is likely that novae 
will be imaged many times, before they erupt, during the course of coming synoptic surveys like 
Pan-Starrs and LSST. We can therefore expect that we will eventually acquire period changes and 
dynamical masses for more classical and recurrent nova ejections. Such data will be invaluable to 
test the methodology proposed in this paper for measuring the ejected masses of all novae with 
well measured distances and UV through IR light curves, and especially to help constrain models 
of SNIa progenitors. 
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5. Conclusions 

We have demonstrated that the assumption of equipartition of energy between radiation and 
kinetic energy in nova ejecta is incorrect. The ratio E^/E^ for novae varies by four orders of 
magnitude. Any nova masses derived under the assumption of equipartition of energy may thus 
be wrong by up to a factor of 10,000. We find that there is no correlation between -E , ra d/-E'km and 
either Y or Z. However, E^d/E^ is well-correlated with the average velocity of ejection v during 
a nova eruption. This leads to the remarkable result that the total ejected mass from a nova is 
directly proportional to, and dependent only on, the total bolometric energy radiated by that nova. 
A simple formula for the ejected mass in a nova explosion is now available if the nova distance and 
bolometric light curve can be measured. 

We thank an anonymous referee for excellent suggestions. 
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Fig. 1. — Ratio of radiated to kinetic energy of a nova outburst obtained from 58 model calculations 
for several parameter combinations, as marked. Different symbols correspond to different accretion 
rates, given in the legend in units of logM (Mq yr _1 ). Red, green and blue symbols correspond 
to core WD temperatures Iwd of 10, 30 and 50 xlO 6 Kelvins. 
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Fig. 2. — The ratio of kinetic to radiated energy in 58 nova models as functions of helium abundance 
Y and metallicity Z. No correlation is seen. 
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Fig. 3. — The ratio of kinetic to radiated energy in 58 nova models as a function of average ejected 
velocity. The clear correlation is discussed in the text, and is key to determining the masses of nova 
ejecta. 
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Fig. 4. — The ejected mass as a function of radiated energy in 58 nova models. The clear correlation 
is discussed in the text; this is our key prediction for determining the masses of nova ejecta. 



